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INTRODUCTION 

The longstanding technological demand for improved 
performance in materials continues unabated. In fact, the 
promise of more than a decade in filamentary composite ma- 
terials has generated intense interest. Such interest is 
well founded, having been generated by the possibility of 
using brittle materials with high modulus, high strength, 
but low density in composites which fail in a non-cata- 
strophic manner. Such fiber reinforced composite materials 
offer substantially improved performance and potentially 
lower costs for aerospace hardware. 

Fulfilling the promise of composite materials, however, 
requires a strong technology base. NASA and AFOSR have re- 
alized that to fully exploit composites in sophisticated 
aerospace structures the technology base must be improved. 
This, in turn, calls for expanding fundamental knowledge 
and the means by which it can be successfully applied in 
design and manufacture. It also calls for expanding the 
body of engineers and scientists competent in these areas. 

As part of their approach to accomplishing this, NASA and 
AFOSR have funded the current composites program at Rens- 
selaer. The purpose of the RPI composites program is to 
develop advanced technology in the areas of physical prop- 
erties, structural concepts and analysis, manufacturing, 
reliability and life prediction. Concommitant goals are 
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to educate engineers to design and use composite materials 
as normal or conventional materials. A multifaceted pro- 
gram has been instituted to achieve these objectives. The 
major elements of the program are: 

1. CAPCOMP (Composite Aircraft Program Compo nent) . 
CAPCOMP is primarily a graduate level project being con- 
ducted in parallel with a composite structures program 
sponsored by NASA and performed by a private, aerospace 
manufacturing contractor. The first component redesign is 
being done in conjunction with the Boeing Commercial Air- 
plane Company. The main spar/rib region on the Boeing 727 
elevator, near its actuator attachment point, has been se- 
lected, with Boeing's advice and the concurrence of NASA/ 
AFOSR, for study in CAPCOflP. The magnitude of the project 
- studying, designing, fabricating and testing a relatively 
small but highly stressed region on the elevator - is both 
consistent with Rensselaer's capabilities and a significant 
challenge. The selection of a portion of a full-scale 
flight hardware structure assures relevance to this pro- 
ject's direction. 

Visits to Boeing were conducted in the Fall of 1978 
by Professor Hoff and several of his students, and the first 
serious design work began shortly thereafter. Two alterna- 
tive designs were pursued to the point of preliminary 
analysis and testing; further progress is reported in Part 


I. 
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2. CAPGLIDE (Composite Aircraft Program Glider ) . 

This undergraduate demonstration project is to design, fab- 
ricate and test a foot-launched, ultralight glider using 
composite structures. A flight vehicle was selected to 
maximize student interest and to provide the students with 
a broad-based engineering experience. For those students 
continuing with graduate work at RPI , CAPGLIDE is intended 
to provide natural progression to CAPCOMP. The progress 
on the CAPGLIDE project to date has been very good. Seven 
professors and approximately 30 students were actively en- 
gaged in the project during the beginning of this period; 
that is, at the end of the spring semester. A description 
of the status of the work performed under CAPGLIDE to the 
end of this reporting period, .is given in Part II. 

COMP AD ( Comp uter Aided Design) . A major thrust 
of the composites program is to develop effective and effi- 
cient tools for the analysis and design of composite struc- 
tures. Rensselaer and NASA Langley have jointly implemented 
the use of the SPAR code on minicomputers, and the work at 
Rensselaer has made "/xrtual memory" available to those 
using SPAR. Recently attention has turned to preprocessor 
developments; details are reported in Pat . CII. 

4. Composites Fabrication and Test Facility . Struc- 
tural design engineers, educated only by course work and 
design projects limited to paper, often fail to sense or 
appreciate problems involved in fabrication. The actual 
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fabrication and testing of composite structural components 
provides this training and the final validation for the de- 
signs in our CAP projects. RPI's Composites Fabrication 
and Test Facility is located in the laboratory and high bay 
areas of the Jonsson Engineering Center. Equipment is 
available for comprc. sion molding parts as large as 19” x 
19" and vacuum bag.jing parts of much larger size. Panels 
approximately 4' x- 20' have been made by vacuum bagging. 
Specifics are given in Part II, CAPGLIDE. A pressure vessel 
for small parts and spars has been designed and built. A 
second one capable of higher pressures is also planned; more 
complete details are reported in Part I, CAPCOMP. 

Our capabilities and efficiency in fabrication and 
testing were increased substantially during the reporting 
period by the addition of a Manager and Master Technician 
for the Composites Laboratory. This position is filled by 
Volker Paedelt whose background includes almost twenty years 
at the Vereinigte Deutscne Metallwerke, Frankfurt, Germany; 
Simmons Fastener, Menands, N. Y. and the Physics Department 
at Rensselaer Polytechnic Institute 

5. INSURE ( Innovative and Supporting Research) . The 
criteria for selection of research projects to be conducted 
under this program are (a) that they must anticipate criti- 
cal problem areas which may occur in the CAP or NASA/AFOSR 
programs or (b) that solutions to existing problems are not 
yet satisfactorily in hand. During the reporting period 
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five programs were funded. Results from the ongoing projects 
are reported in Part IV. 

Curriculum Revisions . The goal of educating engi- 
neers to think of composites as normal or conventional mate- 
rials has required changes in curriculum. Since the initi- 
ation of this program, almost all Rensselaer engineers take 
introductory courses which incorporate the concepts of ani- 
sotropy and composite materials. In addition, six special- 
ized courses in composites have been offered during the past 
three years to develop those special skills required of 
students involved in the composites program. A new course 
was introduced in the Fall of '78 semester on composite 
design and analysis using programmable hand calculators, a 
central mini and full frame computers. A new graduate level 
advanced topics course with the title "Advanced Finite Ele- 
ments" was offered for the first time in September 1979. 

The additions of the SAP and SPAR computer codes and 
the growing availability of interactive computer graphics 
under our COMPAD program element have reached the point 
where our engineering students are using these facilities 
as everyday working tools for design, analysis and visuali- 
zation purposes. We have thus achieved one of the principal 
goals of the curriculum development activities. 

7 . Technical Interchange . 

a) Dr. Kiyoshi Kenmochi returned to the Industrial Products 
Research Institute of Japan's Industrial Science & Technology 
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Agency in August. He was most helpful during his appoint- 
ment as Research Associate under the program, having set up 
our photoelastic stress analysis facilities and assisting 
us in conducting the mechanical attachment research reported 
in our last report (July, 1979) and herein, 
b) Technical Meetings: Technical meetings, on- and off- 

campus, provide important opportunities for interchange of 
technical information. Because of the large number of com- 
posites meetings, a central catalog with all upcoming meet- 
ings is being maintained and distributed periodically. In 
this way we help assure that a Rensselaer staff member will 
participate in important meetings. The calendar for this 
reporting period is shown in Table I. Meetings attended by 
RPI composites program faculty/staff during the reporting 
period are shown in Table II. Some meetings particularly 
relevant to composites, held on-campus with off-campus 
speakers, are listed in Table III. A list of composite-re- 
lated visits to relevant organizations by RPI faculty/staff/ 
students, with the purpose of each visit outlined, is pre- 
sented in Table IV. Regular readers of these progress re- 
ports will note the higher level of such activity in the 
months from September to April - partially because full aca- 
demic year programs are in session - than in those reports 
covering periods including the summer. 

In summary, the NASA/AFOSR Composites Aircraft Program 
is a multi-faceted program whereby aeronautical, mechanical 
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1979 

4/4-6 

4/16-27 

5/7-10 

5/9-11 

5/14-17 

5/21-25 

5/21-25 

5/22-23 

5/23-25 

6/14-16 

6/18-20 

6/25-29 

6/25-29 

7/11-13 

8/6-10 

8 / 20-22 

9/3-7 

9/10-12 

9/10-13 

9/11-14 


TABLE I 

CALENDAR OF COMPOS ITETj -RELATED MEETINGS 
for the period April '79 to September '79 


20th Structures, Structural Dynamics 6 Materials 
Conf., St. Louis, Mo. Sponsored by AIAA/ASME/ASCE . 

Very Large Vehicle Conference, Arlington, Va. 

Sponsored by AIAA. 

Design Engineering Conference & Show, Chicago, 111. 
Sponsored by ASHE. 

Materials Conference, Washington, D. C. Sponsored 
by ASME. 

Conference on Space Manufacturing Facilities, Prince- 
ton, N. J. Sponsored by Princeton and AIAA. 

Fracture of Composites, San Francisco, Ca. Sponsored 
by Society for Experimental Stress Analysis. 

nigh Modulus Fibers & Their Composites, San Fran- 
cisco, Ca. Sponsored by ASTM. 

Symposium on the Fatigue of Composite Materials, 

San Francisco, Ca. Sponsored by ASTM. 

Advanced Composites: Design and Applications, 

Gaithersburg, Md. Sponsored by National Bureau of 
Standards, NASA and NADC. 

Reinforced Plastics/ Composites, Chicago, 111. Spon- 
sored by Plastics Seminars. 

Applied Mechanics, Fluids Eng. and Bioeng. Conference, 
Niagara Falls, N. Y. Sponsored by ASME and CSME. 

Pressure Vessels & Piping Conference, San Francisco, 
Ca. Sponsored by ASME. 

Conference on Carbon, University Park, Pa. 

Lighter-than-Air Systems Technology Conference, Palo 
Alto, Ca. Sponsored by AIAA. 

Sixth Annual Conference on Computer Graphics & Inter- 
active Techniques, Chicago, 111. 

Aircraft Systems & Technical Meeting, New York, N. Y. 
Sponsored by AIAA. 

Practical Applications and Design of Fibre Composites, 
University of Cambridge, England. 

Design Engineering Technical Conference, St. Louis, 

Mo. Sponsored by ASME. 

International Conference on Acoustic Emission, Ana- 
heim, Ca. Sponsored by ASNT. 

International Forum: Composites, Assemblage & Adhe- 

sion, Lyon, France. 
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TABLE I X 

COMPOSITES-RELATED TECHNICAL MEETIiN i -iS ATTENDED OFF-CAMPUS 
for the period April '79 to September '79 


1979 

5/21-26 ASTM Symposium on Composite Materials, Visit EPRI 
and Lawrence Livermore Lab. (Prof. Stoloff ) . 

5/21-22 American Helicopter Society "'orura (Prof. Locwy) , 
Washington, D. C. 

5/21-24 ASTM Fatigue Conference on Composites (Prof. Krempl) , 
San Francisco, California. 

6/7-8 Conference at National Bureau of Standards on Tech- 
nical Cooperative Center for Composites (Prof. Le- 
Maistre) , Washington, D. C. 

7/11-13 Acoustical Society of America Meeting, (Prof. Das 
and C. Lanzl) , Boston, Massachusetts. 


TABLE III 

COMPOSITES-RELATED MEETINGS /TALKS HELD AT RPI 
(April '79 - September '79) 


Topic 

Finite Element Analysis of 
the Post Buckling Behavior 
of Curved, Laminated, Com- 
posite Panels 

Structural Design for 
Longivity: A Professor's 

View 

Stress Analysis Around 
Macroscopic Cracks Under 
Creep Conditions 

The Effect of Statistical 
Correlation on Weakest 
Link Systems 

Aromatic Polyamides: De- 

velopment and Uses 


New Methods in High Tem- 
perature Structural De- 
sign 

World Trends in Computer- 
Aided Manufacturing 


Date Speaker (s) 

4/9/79 N. J. Kudva, Virgin- 

ia Polytechnic Inst. 


4/24/79 James W. Mar, Massa- 

chusetts Institute 
of Technology 

4/27/79 Hermann Reidel, Max 

Planck Institute for 
Iron Study 

5/3/79 R. C. Garson, Univer- 

sity of Pittsburgh 


10/10/79 Paul W. Morgan, for- 

merly Research Asso- 
ciate, E. I. duPont 
de Nemours and Co. 

10/11/79 Donald Griffin, Wes- 

tinghouse Corporation 


10/25/79 M. Eugene Merchant, 
Director of Research 
Planning, Cincinnati 
Milacron, Inc. 


TABLE IV 


COMPOEITE-RELATED VISITS TO RELEVANT ORGANI !:ATI0NS 
by RPI Faculty/Staff /Students 


Visited 

Date 


by Prof (s) 

Purpose 

AFOSR; 

5/10/79 

R 

. Loewy 

To discuss NDE , 

Dr. L. Kravitz, 
Director 




Fatigue and 
Failure Mechan- 
isms 

W.P.A.F.B. : 

7/30/79 

J 

. Diefendorf 

To discuss the 

Mr. G. Peterson 
and personnel 
of APML and 
AFFDL 


L 

G 

R 

S 

. Feeser 
. Helwig 
. Loewy 
. Sternstein 

NASA/AFOSR Pro- 
gram and USAFML 
amd USAFFDL pro 
grams in compos 
ites 


and materials engineers must interact to achieve its goals 
"Hard-nosed" engineering of composite aircraft structures 
is balanced against research aimed at solving present and 
future problems. In the following sections, detailed de- 
scriptions of the CAPCOMP, CAPGLIDE, COMPAD and INSURE pro 
grams are presented. 


PART I 


CAPCOMP (Composite Aircraft Program Component) 
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CAPCOMP (Composite Aircraft Program Component) 

(N. Hoff, R. Loewy) 

CAPCOMP is a program to design flight critical struc- 
tures to take the maximum advantage of composite materials. 
By combining the efforts of experienced faculty with bright 
and well trained but inexperienced graduate students in an 
environment relatively free of traditional design and manu- 
facturing processes, we hope to devise new and hopefully 
useful design concepts. 

There is sufficient information available todav to 
prove that many structural elements can be built lighter of 
advanced composites than of metals. But if such elements 
have to be joined by any other method than adhesive bonding, 
difficulties and uncertainties arise which can be eliminated 
only through conservative designs with their attendant pen- 
alties in weight or by extensive and expensive programs of 
"cut and try". This stands as one important impediment ta 
full adoption of composites by the aerospace industry. 

On the basis of these considerations Rensselaer Poly- 
technic Institute began, as the first task aimed at new 
structural concepts, the design using composites of a joint 
used in an airplane elevator. To make the design realistic, 
an existing metal airframe component was chosen for redesign 
in composites. The existing design chosen was that of the 
Boeing 727 elevator actuator attachment. We conceive of 
this work as carrying forward a Structures Demonstration 
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Program using the joint of the 727 elevator, paralleling 
that of NASA and its aerospace engineering contractor, the 
Boeing Commercial Airplane Company. Our design, fabrication 
and test effort will emphasize design ideas specifically 
suited to advanced composite construction for the purpose 
of minimizing the weight of the structure, but on a scale 
consistent with the university context and funding level. 

The staff of RPI is very grateful to the Boeing Company and 
its engineers for their wholehearted support of the work 
at RPI. 

During the reporting period, the two different designs 
reported previously (July, 1979) as suitable for replacing 
the largely metal attachiaent produced by Boeing were pur- 
sued. The first, using quasi-isotropic graphite-epoxy 
laminates (the R. W. Berg design) , was chosen for initial 
fabrication and second-round stress analysis. The second 
(C. Muser) design, in which a deliberate attempt was made 
to use uniaxial graphite-epoxy tape to as great an extent 
as possible, was suspended in favor of new design analyses 
and test planning with the objective of optimizing the 
strength of the movable, heavily loaded bolted joint attach- 
ment point. 

The research aimed at maximizing rhe load carrying 
ability of the attachment of the rib flanges to the elevator 
skin through mechanical (i.e., pin-type) fasteners also con- 
tinued. The literature and theory of empty and loaded 
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circular holes in composite plates was reviewed as part of 
the process of determining the most efficient arrangement 
of fibers around the hole. 

These efforts have been jointly directed by Dr. Nicholas 
J. Hoff, in part-time arrangement, and by Dr. Robert G. 

Loewy. They are described further in the following para- 
graphs . 

1. The Elevator and Its Attachment 

The conventional aluminum alloy elevator of the Boeing 
727 is shown in the upper half of Figure 1. The lower half 
of the figure is the new version of the elevator redesigned 
by Boeing in graphite epoxy; it is evident from the pictures 
that the latter is composed of fewer parts than the former. 
However, the actuator fitting of the new design is still 
manufactured of aluminum alloy. This fitting is shown in 
Figure 2. The fitting is attached to outboard and inboard 
portions of a new graphite-epoxy spar and to a graphite- 
epoxy nomex-honeycomb rib as indicated in Figure 3. 

The attachment was designed by Boeing to carry loads 
up to 19,000 lbs. The direction of the load varies as the 
elevator rotates over an angle of 28 degrees from the full- 
down to the full-up position. 

2. Berg's Design 


Berg's design is shown in Figures 4 and 5. The first 
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BOEING ELEVATOR ASSEMBLY 


Figure 1 









PROPOSED DESIGN (Berg) 


ACTUATOR LUG TO RIB TRANSITION - TOP VIEW 


Figure 4 





Figure 
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of these figures compares the Boeing composite design with 
Berg's design; both are shown. Figure 5 shows how the spar 
and new actuator rib would be assembled and the build-up of 
thickness in the lug areas anticipated, as needed, to provide 
the required bearing strength. Note that the edges of some 
of the layers are bent 90 degrees to form flanges to which 
the upper and lower coverplates of the elevator can be attach- 
ed. Attacliment would be by means of titanium Hi-Loc fasteners. 
The right-hand and left-hand graphite -epoxy webs are stabiliz- 
ed by a layer of nomex-honeycomb between them. 

The finite-element analysis of the actuator rib v;eb 
(including the areas of the lugs with the bolt holes) re- 
ported for the last period has been refined. In the first 
computation, the wall thickness of the actuator rib web was 
taken as constant to examine the stress distribution in the 
area of the holes and the way these stresses transfer loads 
to the flanges of the actuator rib. The mathematical model 
contained 482 nodes and a total of 800 quadrilateral and 
triangular elements (see Figure 6) , introduced forces simu- 
lating bolt loads from the actuator and reaction links (see 
Figure 2) and a distributed shear flow along the upper and 
lower edges of the actuator rib to represent the reacting 
torque resulting from aerodynamic forces on the elevator. 

In the current reporting period the analysis was re- 
peated for a variable thickness web, using the same mathe- 
matical models. Several iterations v;ere made examining 
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stress levels and varying thickness distributions until a 
relatively uniform stress level design resulted. The vari> 
able thickness, finite element model is shown in Figure 7. 
Graphic comparisons of the thickness distribution shown in 
the last report with that resulting fr<m the analysis de- 
scribed h(^re can be made by examining Figures 8a and 8b. 
Similarly, comparison of the stress distribution as last 
reported tor uniform thickness with that calculated as de- 
scribed above can be made by examining Figures 9a and 9b. 

These finite element analyses results revealed two 
possible unconservative features of Berg'., initial design, 
which may otherwise have gone undetected. First, the plate 
thickness around the attachment bolt holes, based on the 
Advanced Composite Design Guide, v/as not large enough to 
insure the desired stress levels. The more detailed model- 
ing described here was necessary to accurately assess these 
stresses. 

Second, there is a large stress concentration at the 
corner just aft of the upper bolt hole. The Berg thickness 
distribution shown in Figure 7 had not yet taken account of 
that, even after the first finite element analysis redesign. 

In preparation for fabricating the first full elevator 
actuator attachment rib to be undertaken, step-by-step plan- 
ning of the fabrication method was begun. This planning is 
to include the tooling, the (quasi-isotropic) composite ma- 
terial lay-up patterns and stacking sequence, the curing 
cycles and estimates of the kind of distortion and residual 
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stresses that would be expected to result using the proposed 
procedure . 

In the first stages of this planning, it oecame clear 
that ease of operation and other efficiencies argue strongly 
for obtaining a small autoclave, just large enough to make 
the kinds of parts that we are considering under CAPCOMP; 
i.e., between one to two feet in diameter by about five feet 
long internally and capable of withstanding pressures in the 
75 to 100 psi range. In the time that it takes to bring this 
new facility on line, experiments with the proposed flange- 
web lay-ups will be undertaken using vacuuia bagging tech- 
niques, to investigate distortion problems. 

3. Optimizing Fiber Orientations in The 
Vicinity of Heavily Loaded Joints 

It has often been noted that a basic advantage of fiber - 
reinforced materials is the possibility of "tailoring" the 
material to transmit forces most efficiently. Since the 
number of possible material choices and fiber directions is 
very large, suitable methods of analysis are being developed 
to advance the methodology for minimizing weight in the de- 
sign of attachments or joints. We have made some strides in 
theoretical analysis and intend to build and test some spe- 
cially designed specimens. Bolted attachments involving 
only a single bolt are being looked at first. 

Our theoretical analysis of single bolted connections 
will consider, at first, not-yet-optimized material arrange- 
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mer.ts near the joint consisting of ordinary orthotropic 
graphite epoxy materials and plates with empty holes. This 
is being done to check the method and to compare our efforts 
to obtain results and the results themselves with those de- 
veloped by others. The next case considered will be the 
loaded hole as found in connections such as the attachments 
of moving control surfaces. The manufacturing problem will 
be kept in mind, but the mechanical efficiency has priority, 
at least for these first stages. 

As reported in previous submissions, the Muser design 
of the actuator rib of the B-727 elevator was based on the 
idea of using fibers running largely tangent to the circum- 
ference of the hole through which the forces from the actuator 
rod would be transmitted into the elevator structure. It is 
hoped that our theoretical work will lead to a new design 
for laminae in the region of the connector. We anticipate a 
lay-up pattern with fibers tangential to the circumference 
near the hole. From experience and from published reports 
(Ref. 12, 19, 20, 21 in Table V) we believe that this arrange- 
ment should result in lower stress concentrations. Our work 
plan has the objective of confirming this theoretically and 
experimentally. 

It should be noted that the method of calculation being 
developed will be suitable only to determine stress concen- 
tration factors in the linearly elastic range. Since bearing 
failure occurs outside this range, calculations valid beyond 
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TABLE V 

LIST OF REFERENCES PERTAINING TO BEARING STRENGTH OF COMPOSITE S 


1. Frocht, M. M. and H. N. Hill, "Stress Concentration Factors 
around a Central Circular Hole in a Plate Loaded through a 
Pin in the Hole", J. Applied Mechanics, March 1940, pp. A5- 
A9. 

2. Dalas, Richard N. , "Mechanical Joints in Structural Compos- 
ites", Advances in Structural Composites, Vol. 12, Society 
of Aerospace Material and Process Engineers, 12th National 
Symposium and Exhibit, October 1967. 

3. Dalas, Richard N., "Methods of Joining Advanced Fibrous 
Composites", Composite Materials; Testing and Design, ASTM 
STP 460, 1969, pp. 381-392. 

4. Schjelderup, H. C. and B. H. Jones, "Practical Influence of 
Fibrous Reinforced Composites in Aircraft Structural Design", 
ASTM STP 460, 1969, pp. 285-306. 

5. Waszczak, J. P. and T. A. Cruse, "Failure Mode and Strength 
Prediction of Anisotropic Bolt Bearing Specimens", J. Com- 
posite Materials, Vol. 5, 1971. 

6. Padawer, G. E., "The Strength of Bolted Connections in 
Graphite/Epoxy Composites Reinforced by Colaminated Boron 
Film", ASTM STP 497, 1972, pp. 396-414. 

7. Zurko, M. J., "Preliminary Investigation for Developing Op- 
timum Fasteners for Joining Composite Materials", Naval Air 
Development Center, AD908930, February 1973. 

8. Greszczuk, L. B., "Compressive Strength and Failure Modes 
of Unidirectional Composites", ASTM STP 521, 1973, pp. 192- 
217. 

9. Kim, R. Y., "Effect of Temperature and Moisture on Pin-Bear- 
ing Strength of Composite Laminates", J. Composite Materials, 
Vol. 10, April 1976, pp. 149-155. 

10. Hart-Smith, L. J., "Bolted Joints in Graphite-Epoxy Compos- 
ites", NASA CR-144899, June 1976. 

11. Quinn, W. J. and F. L. Matthews, "The Effect of Stacking 
Sequence on the Pin-Bearing Strength in Glass Fibre Rein- 
forced Plastic*', J. Composite Materials, Vol. 11, April 
1977, 139-145. 
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TABLE V continued 


12. de Jong, T. , "Stresses around a Pin-loaded Hole in an 
Elastically Orthotropic or Isotropic Plate", J. Compos- 
ite Materials, Vol. 11, July 1977, pp. 313-331. 

13. Kubo, J. T. , "Joints", Advanced Composites Design Guide, 
Vol. II, Wright-Patterson Air Force Base, 1977. 

14. Collings, T. A., "The Strength of Bolted Joints in 
Multi-Directional CPRP Laminates", Her Majesty's Sta- 
tionary Office, London, 1977. 

15. Wilkins, D. J. , "Environmental Sensitivity Tests of 
Graphite-Epoxy Bolt Bearing Properties", ASTM STP 617, 
1977, pp. 497-513. 

16. Eshwar, V. A., "Analysis of Clearance Fit Pin Joints", 

J. Mechanical Science, Vol. 20, Pergamon Press, 1978, 
pp. 477-484. 

17. Fukuda, H. , "Fastening of FRP", Unpublished Report, 
Aerospace Research Laboratory, Tokyo University, 1979. 

18. Collings, T. A., "The Use of Bolted Connections as a 
Means of Joining Carbon Fibre Reinforced Plastics", De- 
sign with Fibre Reinforced Materials, Institute of 
Mechanical Engineers, Conference Publication, 1977-1979, 
pp. 29-32. 

19. Puppo, A. and J. Haener, "Application of Micromechanics 
to Joints and Cutouts", U. S. Army Aviation Material 
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the linear region would probably be necessary, and suitable 
failure criteria would have to be established before there 
can be hope for good agreement between analysis and test re- 
sults. We expect this sort of work to wait till late spring 
because the construction of the test specimens, their test- 
ing and the linear analytical work will require all the time 
available until then. 

4. Supporting Development of Mechanical Joints 

The work of graduate student Wonsub Kim and, until 
August, of Research Associate Kiyoshi Kenmochi, in this pe- 
riod, concentrated on the interpretation of bearing failure 
test and photoelastic results in pin-loaded composite tension 
specimens. This followed from our conclusion that designing 
a mechanical joint to fail in bearing makes the best use of 
strength. Reducing fastener diameter increases a given 
joint's tensile strength and to a first approximation leaves 
shear strength unaffected, but bearing stress is ultimately 
limiting in this process. For this reason, we have concen- 
trated on bearing strength and bearing failure mechanisms in 
this research. 

To place the results we obtained in perspective, a lit- 
erature survey was made of experimental and theoretical work 
published in this area. This survey has shown that v/ithin 
the substantial literature of composite materials, relatively 
few research results connected with mechanical joining prob- 
lems have been published. Furthermore, among those reported 
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works which do exist many are inconclusive and others are 
not mutually consistent. A total of 21 references which 
seem pertinent are listed in Table V on pages 30-31 and a sum 
mary of the most important results, from the point of view 
of this study, follows. 

Summary of the Literature Survey 

A. General 

1) There are relatively few publications 
available dealing with bearing failure 
in composite materials. 

2) Among such publications as do exist the 
work reported is roughly 90% experimen- 
tal, 10% theoretical. 

B. Theoretical VJork 

1) The most pertinent work appears to be 

* 

that of T. de Jong (12) and V. A. 

Eshwar (16) . Both are two dimensional 
theories and assume infinitely rigid 
pins. Neither includes the effect of 
friction, comparisons with experiment 
nor satisfies boundary conditions ex- 
actly. 

2) T. de Jong shov;s marked effects of ply 
orientations on stress distributions, 
but provides no estimates of the in- 
accuracies introduced by failure to 
meet displacement boundary conditions 
on the hole. 


Numbers in parentheses in this section correspond to refer- 
ence numbers in Table V, pages 30-31. 


34 


3) V. S. Eshwar's work is limited to iso- 
tropic materials but, by using an inverse 
procedure specifying displacements and 
then solving for the corresponding loads 
on the hole, shows that the boundary con- 
ditions can be met to within satisfactory 
accuracy. Among his results are the 
following : 

a) The area of contact between pin and 
plate has an asymptotic relation to 
load limited by Poisson's ratio, 
for example, to 166® (compared to 
180®) for V = .3. 

b) Non-linear behavior of maximum 
stresses with load are exhibited 
(concave downward) as a result of 
increasing contact area with load. 

c) The distribution of radial stress, 
however, changes only slightly 
with load and is well approximated 


by a 

r 


0 

r 

max 


1 - 

'0 

2 


0 



where 


0^ is the angular position of the 
edge of contact measured from the 
axis of symmetry. 

d) The effect of the third direction 
is likely to be great near the hole 
and suggests caution regarding use 
of 2-D theory in formulating fail- 
ure criteria. 


C. 


Photoelastic Studies 

Stress concentration factors determined by Frocht 
and Hill (1) were found to vary as follows: 


1) If the ratio of the diamett.r of the hole 
to the width of the plate, D/t increases, 
the stress concentration factor decreases 

2) If the clearance between pin and hole in- 
creases, the stress concentration factor 
decreases . 

3) If the number of pins increases, the 
stress concentration factor decreases. 

4) The effect of lubrication between pin 
and hole is measurable but very small 
(thus validating the assumption of no 
friction between pin and hole in mathe- 
matical analyses) . 

Experimental Work 

1) All referenced works investigate the 
effects of one or more of the follow- 
ing six factors. 

a) Fiber orientation 

b) Stacking sequence 

c) Ratio of hole diameter to plate 
thickness (D/t) 

d) Clamping force 

e) Environmental condition 

f) Local reinforcement of the hole 

2) All tests used virtually the same method, 
i.e., apply load to the pin until failure 
occurs. 

3) Of ten references reporting the effect of 
fiber orientation, four (3, 8, 15, 18) re 
ported the effect to be great, five (2, 9 
10, 14, 17) pronounced the effect small 
and one (13) was noncommital. 
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4) Two (9, 13) of three investigations reported 
environmental factors as important influ- 
ences, but one (15) found them of no signifi- 
cance. 

5) All three of those reporting stacking se- 
quence effects (11, 14, 18) concluded it is 
relatively large (But two of these were 
written by the same author.). 

6) Four (10, 13, 14, 18) out of six papers re- 
ported clamping force as having a large in- 
fluence on bearing strength. The remaining 
two authors (3, 7) pointed out that the 
significant fact is the existence of clamp- 
ing force, not the amount of such force. 

7) Five investigators tested the effect of 
D/t; one (14) reported it as big, two (15, 

17) as small, and two (2, 13) found their 
tests inconclusive. 

8) Three authors (2, 6, 10) investigated the 
influence of local reinforcing around the 
hole using colaminated third materials, 

and they all reported the effect to be favor- 
able and very large. 

9) Two reports (5, 13) mention the effect of 
clearance between pin and hole, but it was 
only calculated using a Finite Element 
Method . 

The results of RPI's testing is presented in this sec- 
tion in some detail and is also related to the corresponding 
data found in the literature search. Table VI summarizes 
RPI bearing test results to date. Analysis of the photo- 
elastic test results hasn't been completed. Some initial 
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TABLE VI 

RESULTS OF RPI BEARING TESTS 


* 

Specimen 

Designation 

Specimen 

Thickness 

(inch) 

Load 

Maximum 

(lbs) 

Bearing 

Strength 

(psi) 

strength Ratio 
to that of 
Specimen E 

** 

Type A 


1,240^^^ 



-1 

0.085 

58,353 

0.82 

-2 

0.083 

1,500 

73,764 

0.94 

Type B 





-1 

0.094 

1,610 

68,510 

0.97 

-2 

0.094 

1,725 

74,902 

0.95 

Type C 





-1 

0.095 

1,445 

60,842 

0.86 

-2 

0.095 

1,625 

69,817 

0.89 

-3 

0.095 

1,120 

47,158 

0.60^^^ 

Type D 





-1 

0.065 

1,400 

65,882 

0.93 

-2 

0.086 

1,310^^^ 

62,173 

0.79 

Type E 





-1 

0.083 

1,470 

70,843 

1.00 

-2 

0.083 

1,600 

78,682 

1.00 

Photoelastic 

0.089 

2,750 

61,798 

0.72 

II II 

0.095 

3,350 

70,526 

0.82 

II If 

0.089 

3,100 

69,663 

0.81 

II II 

0.085 

3,650 

85,882 

1.00 

★ 

Specimen Types A-E have holes 
in the photoelastic specimens 

with h inch 
is h inch. 

diameter . 

The hole diameter 

The specimen 
Type A 
Type B 
Type C 

designations are 
{0°/90°}8s 
{±45®/0®}6s 
{±45°/90®}6s 

interpreted as follows ; 

Type D {±45V0®/90“}4s 

Type E {±45°/±22.5V±67.5®/0°/90°}2s 

^^^Hole drilled 

with high speed 

steel drill ? 

oit . 



( 2 ) 

Found to be premature failure. 
(3) 

Ratioed to Case E-2 . 


t 
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conclusions/ however, can be drawn? namely, 

A. The material behaves as if almost perfectly elastic 
until just in the vicinity of the maximum load. 

Through both the unloading and loading process the 
charting pen follows exactly the same curve of load 
vs displacement. Even after, the unloading and load- 

ing process, the photoelastic fringes were identical " 

to those observed before this process. The behavior 
close to the failure load could not be checked by 
these means. 

B. Stress distributions are quite different from one 
fiber orientation type to another (see Figure 10) . 

On the other hand, the values in Table VI do not 
always confirm the photoelastic insights. "Scatter" 
between the duplicate specimen pairs is large enough 
to almost obscure differences, but specimen Type E, 
which is closest to isotropy, showed the highest 
strength? Type B looks next best? and Types A, C and 
D seem roughly the same and least good. Since hole 
diameter to specimen width varied slightly between 
the photoelastic and other specimens, it may be best 
to consider bearing strength averages based only on 
non-photoelastic tests. Then the strength of lay-up 
Type E is around 75 Ksi. The lowest average would 
be that of Type C at 65 Ksi. This way of looking at 
the data makes the maximum variation of bearing 
strength less than 12%. Hence, contrary to wliat 
would be expected from photoelastio results, such as 
shown in Figure 10, the pin-bearing strength does 
not change greatly from one fiber orientation to 
another. Two exceptions to this statement are speci- 
mens which failed prematurely. 
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Tables VII and VIII relate the RPI results to those in 
the literature. Table IX represents an attempt to draw con- 
clusions as to fiber orientations which constitute good de- 
sign for bearing strength by making comparisions selectively 
among the results reported in the literature and those of 
the RPI tests reported here. Since materials varied from 
test to test, bearing failure loads are normalized to par- 
ticular ply orientations. While hardly conclusive, the fol- 
lowing tentative conclusions can be drawn: 

A. Combinations including 0® and ±45* plies have 
superior bearing strength. (Comparisons 2, 3 
and 4 in Table V) 

B) Having 90® plies provides questionable benefits. 
(Comparison 2 and 4) 

C) Clamping does not alter the relative strength 
between 0/90 and 0/±45 plies. (Comparison 2) 

D) Clamping does seem to alter the relative 
strength between 0/±45 and 90/±45 plies. (Com- 
parison 4) 

E) Changing fiber material may alter the relative 
strength between 0/90 and ±45 plies. (Compari- 
son 1, assuming conclusion C) holds also for 
those ply comparisons) 

An interesting conclusion from our testing is that the 
bearing strength does not seem to be a strong function of the 
ratio of the diameter of the hole to the thickness of the 
plate. This is contrary to results presented in Reference 
14. This D/t ratio for photoelastic specimens was roughly 



THESE Nli-IBERS REFER TO THE LIST OF PAPERS IN T^B! E V 
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TABLE IX 

NORMALIZED BEARING STRENGTH COMPARISONS AMONG REFERENCE DATA 

AND FROM 

RPI WORK TO DETERMINE THE EFFECT OF FIBER ORIENTATIONS 


Comparison 

Number 


Fiber 

Orientation 


1. 

Ref . No . — 

* 

14 

17 



Mat'l. 

G/E 

Glass/E 


{0/90}s 


1.0 

1.0 


{±45}s 


1.04 

1.2 


2. 

Ref . No . — 

* 

14 

RPI 



Mat • 1 . 

G/E 

G/E 


{0/90}s 


1.0 

1.0 


{0/±45}s 


1.12 

1.17 


3. 

Ref. No. — 

* 

14 

** 

6 



Mat ' 1 . 

G/E 

G/E 


{±/45}s 


1.0 

1.0 


{0/±45}s 


1.08 

1.7 


4. 

Ref. No. — 

9 

15 

RPI 


Mat ‘ 1 . 

G/E 

G/E 

G/E 

{0/±45}s 


1.0 

1.0 

1.0 

{90/±45}s 


0.88 

0.51 

0.89 

{0/90/±45}s 


1.06 

— 

0.97 

"k k it 

3.2 Ksi clamping; 200 lb 

*** 

clamping; " 

5 lb- 

in on a screw" 


clamping. 
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5.6 and for the bearing coupons was 2.8, yet the correspond- 
ing bearing strengths were quite close. With such ratios, 
according to Reference 14, the bearing strength of the photo- 
alastic specimen should be lower than the bearing test coupon 
by at least 50%, but this series of tests at RPI shows that 
the photoelastic specimens have strengths which are even 
slightly higher. 


5. Failure Mechanisms 

Some of the specimens tested to failure in our program 
were cut along the center line passing through the two holes 
so as to allow microscopic investigation through the speci- 
men thickness in the area of the failure. 

After very careful polishing, the cut edges were in- 
spected both on the failed area of the hole as well as on 
the unfailed side. Microscope magnification factors of 25 
and 50 were used. With the magnification factor of 50, the 
cracks on the damaged area as well as each layer of fiber 
orientation could be seen clearly. Figure 11 shows the 
cracks on failed specimens with several typical fiber orien- 
tations. These pictures were taken with a magnification 
factor of 25. 

It seems clear that the cracks started at an arbitrary 
point on the loaded edge of the hole and ran across the lay- 
ers, even through the 0“ layers, in a direction roughly +45® 
to the laminae. This kind of crack has been observed often 
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in tests of brittle isotropic materials loaded in compres- 
sion. If these cracks are along shear lines as they appear 
to be and as they are in compression tests of isotropic ma- 
terials, then shear in a direction transverse to the plane 
of the plate should be investigated and regarded as an im- 
portant factor in the bearing failure of composite materials. 

6. Special Delamination Studies 

The effect of delamination on the bearing strength of 
a composite material seems to be much more severe than it 
is on other material allowables. In a brief investigation 
of this effect two additional specimens were prepared using 
exactly the same procedures as those used in all the other 
tests. In this case, hov/ever, the holes were purposefully 
drilled with a high speed steel drill bit instead of a car- 
bide tipped drill bit. Every dimension and process of prep- 
aration was identical except the drill bit used. As expected 
these holes were flawed to a greater extent than those drill- 
ed with a carbide tipped bit. The thickness of the plates 
around these holes was increased by up to one thousandth of 
an inch. This indicates that delaminations were induced by 
drilling with the high speed steel bit. These holes were in- 
vestigated further using ultrasonic NDE techniques. 

Figure 12, showing ultrasonic transmission loss images 
of the delamination around holes drilled with a high speed 
steel bit and with a carbide tipped bit, shows the former to 
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be more extensive. Note that the technique does not yet 
allow interpretation of delamination as a function of depth, 
but only of area around the hole in the plane of the lami- 
nae. 

These two specimens were tested using the same process 
as described above. A mishap failed one specimen of lay-up 
Type A ({0*/90®}gg) before useful data could be taken. Spec- 
imen Type C ({±45?/90®}gg) , however, was tested successfully. 
It showed much less strength than the same specimen with 
higher hole quality, as expected. This result is also in- 
cluded in Table II. 

These special tests lead to the following conclusions: 

1) Holes drilled with carbide tipped bits 
have delamination like that caused by 
high speed steel bits, only less of it. 

2) The amount of delamination caused by 
drilling affects the bearing strength 
of a joint substantially. 

3) Preventing or compensating for delami- 
nation in some way can improve the bear- 
ing strength of a joint greatly. The 
large and favorable effect of clamping 
seems to stem from this phenomenon. 


PART II 


CAPGLIDE (Composite Aircraft Program Glider) 
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CAPGLIDE (Composite Aircraft Program Glider) 

(R. J. Diefendorfr H. J. Hagerup^ G. Helwig) 

CAPGLIDE is an undergraduate program to design, build 
and test advanced composite structures. Students will ob- 
tain direct "hands-on" experience in advanced composite 
structures which can serve as a springboard for the more 
sophisticated CAPCOMP projects. In dealing with the design 
of a complete vehicle, the effect of any given change on 
other aspects must be dealt with. In this way the project 
also requires students majoring in aeronautical, mechanical 
and materials engineering to interact in much the same way 
as they do in industry. 

An ultra-light sailplane was selected as the first 
demonstration project because a full scale flight vehicle 
would maximize student interest and would be of relative 
simplicity and lov/ cost to build. A conventional layout 
monoplane with three-axis control resulted in the following 
estimated performance: 

1) Stall speed, 15 knots, 

2) Best glide ratio, 17 and 

3) Minimum sink rate, 2.0 feet per second. 

While the glide ratio of the ultra-light sailplane is simi- 
lar to that of post World War II utility gliders, the more, 
important sink rate is in the range of standard class sail- 
planes . 
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The major activity during the current reporting period 
was the manufacture and assembly of the glider. The fuselage- 
booms were built first. Parallel to this, the wing-skins were 
la id-up on a flat polyethylene-surfaced table and subsequently 
curved to shape in wing molds. This two step method (see 
Figure 13), consisting of (1) semi-curing on a flat table and 
(2) final curing into shape in molds, was very successful. 


Step 1. Plat lay-up on polyethylene plate. 


to vacuum pump 


vacuum bag sealant 



1 


vacuum bag 
foam 




polyethylene plate 




wet laminate 


Step 2. Lay-up in mold. 


vacuum bag 


vacuum bag sealant 


q sealant X ^ / ^foam 





iJIlJ 




vacuum pump 



lastic sheet 
wet laminate 
foam laminate from 
1, laminate 
faces down 


pressed down into mold 
by vacuum 


Wet-Lay-Up Mold Technique 
Figure 13 
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PAr.r 

•' Ot */\i 

Before the cap-strips of the wing spars wore manufac- 
tured, a test spar, which was built earlier, was loaded to 
destruction. A test loading jig was hung from a crane and was 
used to introduce the forces at four points into the canti- 
levered spar (see Figure 14). 


Figure 14 Spar Test 

The cantilever point was achieved by clamping with the con- 
nection from the other wing. The wing-connection was, there- 
fore, tested in the same process. Special guides were neces- 
sary to prevent the I-beam spar from twisting, and these can 
also be seen in Figure 14. In the first attempt, the loading 
jig failed at 50% of the ultimate load (i.e., derign load 
with safety factor of 1.5), and the test spar v/as damaged and 
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had to be repaired. The next test was successful. The spar 
failed in the connection with a safety factor slightly over 
1.5. The nature of the failuie was graphite capstrips split- 
ting under shear forces. The actual wing connection will be 
reinforced to prevent shearing-off of the I-beam flanges. A 
still higher safety factor can* therefore, be expected. A 
3-foot wing section, previously tested for wing-skin buckling, 
was tested again at an elevated temperature of about 125®P. 
Significant stiffness loss was obvious. Accordingly, the 
wing was reinforced with two additional webs from the root 
to the inboard end of the ailerons, which added some extra 
^ weight. 

The tail was built with single-sided, Kevlar-stabilized 
foam sheets. These sheets were then bent to the rib frame- 
work. This technique, again, was simple and effective in 
producing true-to-tolerance airfoil shapes with an extremely 
lightweight structure. 

From the following Figures, 15-18, the construction 
details and materials actually used in fabricating the glider 
can be seen. The actual weight of the parts is also listed. 
Weights of major components sub-assemblies are shown in Table 
X. Figures 19-21 are photographs which show various assem- 
blies of the first CAPGLIDE aircraft, RP-1, in different 
stages of completion; Figure 19 is the tail assembly. Figure 
20 is a semi-span wing substructure and Figure 21 is the 
fuselage. 
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Item Units 


Span 

m(ft) 

11.5 

(37.7) 

Length 

m (ft) 

5.78 

(18.9) 

Height 

m (ft) 

1.28 

(4.2) 

Wing area 

m^(ft‘^) 

12 

(129.0) 

Aspect ratio 


11 


Wing section profile 


FX-63-137 


Vertical stabilizer ~ - 

Area 

m^(ft^) 

0.72 

(7.8) 

Span 

m(ft) 

1.2 

(3.9) 

Profile 


FX-L-111-142 


Horizontal stabilizer ^ ~ 

Area 

m^(ft^) 

1.8 

(19.4) 

Span 

m (ft) 

3.0 

(9.8) 

Profile 


FX-L-111-142 


Performance with 120 kg 

Stall speed 

km/h (mph) 

36 

(22.5) 

Minimum sink rate 

m/s (ft/s) 

0.62 

(2.0) 

at 

km/h (mph) 

40 

(25.0) 

Maximum glide-ratio 


20 


at 

km/h (mph) 

52 

(32.5) 

Empty weight 

kg (lbs) 

52.9 

(116.3) 

Tail^control effectivenesses 

V/V Horizontal 


0.41 


V/V Vertical 


0.057 







Part 

Materials 

Weight (c 

2 Spar capstrips 

Graphite 

1 , 600 

Main web 

Kevlar, foam 

700 

17 Ribs 

Glass, foam 

1,020 

2 Skins 

Kevlar, foam 

CO 

O 

Front web 

Glass, foam 

14 0 

Reinforcement webs 

Glass, foam 

340 

Back web 

Glass, foam 

83 

Root rib 

Graphite, foam, spruce 

750 

Wing connection 

Graphite, plywood, foam 

500 

Aileron tube 

Graphite 

220 

Hinges 

Graphite 

100 

Wing tips 

Kevlar, foam 

900 

Leading edge 

Balsa 

500 


13,433 

2,600 


Total = 16,033 


Resin for bonding 
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Part 

Materials 

Weight (g) 

4 Boom capstrips 

Graphite 

3,000 

4 Wobs 

Kevlar , foam 

3,400 

Stiffeners 

Kevlar, foam 

200 

Fuselage 

Kevlar 

900 

Pilot beam 

Graphite, balsa 

1,300 

Landing skid 

G3.ass, spruce 

1,100 

4 Wing supports 

Spruce 

320 

Push Rods 

Aluminum 

500 

Aileron Tube 

Graphite 

80 

Cranks 

Graphite 

100 

Control stick 

Graphite 

120 

4 Landing struts 

Graphite, balsa 

600 

2 Connections blocks 

Spruce 

140 

Torsion box 

Graphite, foam 

810 

3 Boom connections 

Graphite, balsa 

200 

8 Filler blocks 

Spruce, balsa 

280 

13,050 

Resin for bonding 


2,500 


= 15,550 


Total 
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Figure 18 RP-1 

TAIL-CAPGLIDE 

niNasaLAKii polvtkchnic 

INSTITUTB NOVIMBCn ia7a 


Part 

Materials 

Weight (g) 

Skins 

Kevlar, foam 

3,000 

24 Ribs 

Foam 

144 

Web (horizontal) 

Kevlar, foam 

192 

Web (vertical) 

Kevlar, foam 

96 

2 Capstrips 

Graphite. 

225 

4 Rudder Webs 

Foam 

54 

11 Hinges 

Graphite 

220 

Rudder connection 

Graphite, tube 

150 

Leading edge 

Balsa 

135 

Wing tips 

Balsa 

90 

2 T-Webs 

Kevlar , foam 

128 

Tail connection 

Graphite 

50 


4,484 

600 


Adhesive 


Total 


5,084 
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table X 

CA PGLIDE j RP-1 ACTUAL TOTAL WETnuTC 



(kg) 

(lbs) 

Wings 

32.2 

70.8 

Tail 

5.1 

11.2 

Fuselage 

15.6 

34.3 

Total Weight 

52.9 

116.3 



Tail Assembly 






iry 


Figure 19 
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Figure 21 


Fuselage 




PART III 


COMPAD (Computer Aided Design) 
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COMP AD (Computer Aided Design) 

(L. J. Feeser and M. S. Shephard) 

During this reporting period, emphasis within the com- 
puter aided design portion of the composites project has 
been put on enhancing RPI's finite element capabilities. 

This work is primarily concerned with three activities: 

(1) extension and maintenance of the finite element analysis 
programs available on RPI's computer systems, (2) initial 
development of a highly interactive finite element prepro- 
cessor and (3) linking the finite element analysis program, 
SPAR, to an array processor. 

1. Finite Element Analysis Programs 

Work has begun on the development of a new in-house 
finite element analysis program. This work is being carried 
out primarily as a series of term projects by a few selected 
students in RPI's finite element courses. This new program 
is being developed in a highly modular form so that new cap- 
abilities can be added with a minimum of effort. 

Several other finite element method analysis programs 
are available to us; in fact, enhancement of RPI's finite 
element analysis capability has been accomplished mainly by 
the acquisition of additional finite element analysis codes. 
Table XI lists the finite element programs currently avail- 
able along with a brief description of their capabilities. 
Nevertheless, we anticipate that at some point in the future 
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TABLE XI 

LIST OF FINITE ELEMENT CODES AVAILABLE ON RPI'S COMPUTER SYSTEMS 


Program 

SPAR 

SAP IV 
NONSAP 


ASHSD2 


STRUDL 

ANSYS* 


Computer 

PRIME 

IBM, 

PRIME 

IBM 

IBM 

IBM 

PRIME 


Capability 

Static and dynamic analysis 
of both linear and nonlinear 
structural systems. Especi- 
ally well suited for compos- 
ites. 

Static and dynamic analysis 
of linear three-dimensional 
structural systems. 

Same as SAP IV but with non- 
linear capabilities. Nonlin- 
earities may be due to large 
displacements, large strains 
and nonlinear material behav- 
ior. 

Evaluates the time dependent 
displacements and stresses of 
complex axi symmetric struc- 
tures subjected to any arbi- 
trary static or dynamic load- 
ing or base accelerations. 

Static and dynamic analysis 
of linea’- structural system. 

A wide variety of capabili- 
ties are available in this 
program including: 

a) static and dynamic analy- 
sis of linear and nonlin- 
ear systems (elastic, 
plastic, creep, swelling 
and large deflections 
handled) , 

b) steady-state or transient 
heat transfer analysis 
(conduction, convection 
and radiation) and 

c) some fluid and electrical 
capabilities. 


By contract the use of ANSYS is limited to ''small'' problems 
for course work and student projects. It will not be avail- 
able until Fall '79. 
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we will need to develop some new, then unavailable capahil - 
ites, and it will be far more efficient to add these capabil- 
ites to a general program that we know intimately, than it 
would be to modify one of the other existing programs or to 
write a complete program. An additional enhancement which 
will be discussed in conjunction i^ith the array processor 
work is the creation of a version of SPAR that will run on a 
POP 11/40 minicomputer. 

2 . Preprocessor Development 

The finite element preprocessor currently under develop- 
ment will handle virtually all possible two-dimensional to- 
pologies. Extensive use of graphical input techniques and a 
unified overall design will result in a preprocessor that is 
both efficient and simple to use. 

The development of a finite element model using a general 
preprocessing capability can be thought of as a three stage 
process. The first and most general phase is defining and 
discretizing the problem geometry. The second phase is de- 
fining problem-specific attributes such as material proper- 
ties, boundary conditions and loads. The third phase is 
linking the preprocessor with the particular analysis pro- 
gram selected to solve the model. With this approach, only 
one general program for defining and discretizing problem 
geometry is required, since this process is the same for all 
finite element models. This program can then be linked to 
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one attribute editor for each class of problem to be con- 
sidered. For example r one attribute editor may be used for 
stress analysis problems , another for heat transfer analysis 
and a third for fluid flow calculations. The last stage, 
the link to an analysis routine, is little more than a sim- 
ple reformatting program; it is the only portion of the 
preprocessor program that will have to be repeated for each 
new analysis routine considered. 

A two-dimensional preprocessor has been developed to 
the point that basic problem geometries can be defined and 
element meshes can be generated using a set of mesh gener- 
ating capabilities. The operation of the preprocessor will 
be demonstrated in terms of a simple example. 

The first step is the definition of a background plan- 
ning grid which can be thought of as a piece of electronic 
graph paper. This planning grid is used later to obtain nu- 
merically exact input with interactive input devices by "hit 
testing" on the intersection of grid lines. In addition to 
defining and scaling a basic rectangular grid, interactive 
features allow the user to add and delete additional vertical 
and horizontal grid lines at user-specified levels and to 
move the origin to any grid intersection to produce a final 
planning grid such as the one shown in Figure 22. 

The second step is to define the subregion boundary 
curves to be used by the mesh generators. Curves are defined 
interactively by first selecting the type of curve to be 
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generated and then inputting the appropriate number of "key 
points” (either by use of the light pen, thumb wheels or 
keyboard) along the curve. The curve-generating options 
include single points, lines, circles, circular arcs and 
second and third degree polynomials computed by Lagrangian 
interpolation, A specified number of "key nodes" are then 
placed along the curve. The key nodes are the nodal loca- 
tion used by the mesh generator to define the resulting mesh. 
Since the number and location of the key nodes define the re- 
sulting mesh, options are available to specify both the num- 
ber of key nodes to be placed on a curve and also a weighting 
parameter that forces the nodal spacings to vary from uniform 
along the length of the curve to being biased toward either 
end. Figure 23 shows the interactive input of a circular 
arc by three points along the arc while Figure 24 shows the 
final set of subregion boundary curves for the example prob- 
lem. 

After the subregion boundary curves have been defined, 
the geometry of the object has been determined and the element 
mesh can be generated. The element mesh is generated one sub- 
region c.t a time by employing one of .several available mesh 
generators. 

Currently, three mesh generators have been implemented. 
The generators are computationally efficient and will produce 
well-conditioned grids in highly distorted regions. Each 
generator is based on a mathematical mapping that defines 
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internal node point= in terms of the appropriate boundary 

key nodes by linear blending interpolants in one, two or 

* 

three directions . After the locations of the internal nodes 
have been defined, elements can be specified by appropriately 
connecting the nodes. 

The mesh generators are invoked by first interactively 
indicating the subregion boundary curves that define the 
sides of a subregion, selecting an element type and then in- 
dicating the mesh generator to be used. Considering the left 
hand subregion of the example problem, the user indicates 
that the first three sides are one curve each, while the 
fourth side (the bottom of the subregion) is the concatena- 
tion of two curves. The user then indicates that he wants a 
four-sided subregion generated by the linear blending func- 
tions using four noded quadrilateral elements to define the 
mesh. The resulting mesh is shown in ]‘igure 2 5 where a zoom 
feature has been used to look more closely at that subregion. 
Triangular elements can also be generated as shown in Figure 
26. Figure 27 shows triangular elements generated in a 
three-sided region. 

The upper right hand portion can be generated as a single 
subregion by mapping, via the linear blending functions, from 


* 

Haber, R. B., M. S. Shephard, J. F. Abel, R. H, Gallagher 
and D. P. Greenberg, "A Generalized Graphic Preprocessor 
for Two-Dimensional Finite Element Analysis", Computer 
Graphics a quarterly report of SIGGRAPH-ACM, Vol. 12, No. 
3, August 1978, pp. 323-329. 
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the circle defining the hole to the three curves defining 
the outside of the subregion. Figure 28 shows a four element 
deep mesh, biased towards the hole, generated in this manner. 
Figure 29 shows the completed mesh for the example problem. 

3. Array Processor for SPAR 

As discussed in the last semi-annual progress report, 
preliminary feasibility studies have indicated the potential 
benefits of using an array processor in conjunction with the 
SPAR finite element analysis program. Work is currently 
under way to link SPAR with a Floating Point Systems AP120 
array processor. Since this array processor is operated 
from a POP 11/40 mini-computer, effort has concentrated on 
creating a modified version of SPAR to run on the 11/40. As 
this work is being conducted, areas of code where an array 
processor may be used are being identified and considered as 
candidates for replacement by array processor operations. 



Figure 28 

Two-Sided Mesh 
of Three-Noded 
Elements 


Figure 29 
»csH il Final Mesh 
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PART IV 

INSURE (Innovative and Supporting Research) 


Progress on composites research is reported in the indi- 
vidual write-ups on the succeeding pages in the following 
areas : 

Advanced Analysis Methods for Composite Structures 
Ultrasonic Non-Destructive Testing Developments 
Physical Properties of Epoxy Resins and Composites 
Fatigue in Composite Materials 

Transverse Thermal Expansion of Carbon/Epoxy Composites 
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ADVANCED ANALYSIS METHODS FOR COMPOSITE STRUCTURES 

Senior Investigator; E. J. Brunelle 

During the reporting period three diverse areas of re- 
search were undertaken. The first two areas appear to be 
relatively short term ventures even tnough one of them has 
become diffuse enough to lead to a number of sub-areas that 
are each producing interesting results. The third area is a 
more ambitious undertaking and is considered to be a longer 
term effort with potentially promising and useful results. 
Interestingly enough, research in the first two areas was 
initiated only after assuming that the results should have 
been in the literature; however no trace of them was found, 
and it was, therefore, deemed worthwhile to spend some time 
filling in the theoretical gaps. 

The first area of research delineates the laws governing 
the variation of the properties of an orthotropic lamina, of 
an anisotropic lamina and of a laminate comprised of ortho- 
tropic laminas with changes in the orientation of the coordi- 
nate axes. A nondimensional number B = completely governs 

these rotational properties of an orthotropic lamina and de- 
pending on whether B>1 or B<1, "regular" or "anomalous" rota- 
tional properties are observed. This is not simply an "aca- 
demic" issue since typical values of B for boron/epoxy, 
graphite/epoxy and glass/epoxy are .99, 1.01 and 1.3 respec- 
tively. Besides obtaining detailed rotational information on 
the (elements of the reduced stiffness matrix in the 
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transformed stress-strain relation) as influenced by B, a 
simple inequality has been derived that identifies a new com- 
posite material as being "regular" or "anomalous". This in- 
equality (see Figure 30) should be of particular interest to 
the physical chemists/materials engineers who are involved 
with developing new fibrous composite materials. 

The second area of research traces its origins to a de- 
sire to reduce the number of elastic constants in any given 
problem to a smaller number of "cluster constants" which can 
be used more efficiently in parameter studies of any given 
problem. For example, if a problem involving four elastic 
constants could be reduced to a problem involving only one 
"cluster constant", physical insight into the problem becomes 
much clearer, more general solutions are usually possible. 
Furthermore, what was originally an enormous computational 
effort (Note that for a parametric study of the above cited 
generic problem, sufficient sets and combinations of the four 
constants must be considered.) can be reduced to a one para- 
meter family of solutions. 

"Cluster constants" have been derived by employing af- 
fine transformations on the dependent variables as v/ell as 
the independent variables. This method has been quite suc- 
cessful and has produced some diverse results which follow; 

1) A composite shell problem that was posed, but 
not solved correctly, alraost two dticades ago 
in a classic paper b” Pister, Dong and Taylor 
is now easily solveo 


0.08 0.16 0.24 
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2) The dependence of shear deformation effects 
on boundary conditions has been shovm drama- 
tically by a general solution to a composite 
plate strip problem. The result, which is 
well documented in earlier work for (simpler) 
transversely isotropic beams and plates, is 
that the effects of shear deformation are 
exacerbated by any form of increased bound- 
ary restraint; obviously this is not the 
trend one would prefer to see. 

3) The specially orthotropic plate, which pro- 
vides important bounds for composite plate 
problems (upper bounds for buckling, etc.), 
ca i be reduced to a one parameter family of 
solutions. Master curves can thus be pre- 
sented for any and all static, dynamic and 
instability problems associated with the 
specially orthotropic plate. 

The third area of research seeks a systematic method of 
generating hierachies of (elastic) composite beam, plate and 
shell equations that govern the static, dynamic and stability 
behavior of realistic composite materials structures. In 
this development, a proper accounting is to be made of (i) 
conservative and non-conservative arbitrary initial stress 
states, (ii) the accompanying .’nitial displacement gradient 
field [this has been almost completely neglected in the lit- 
erature of even single layer structures] and (iii) arbitrary 
states of initial [stress-free] displacements so that the 
effects of "imperfection sensitivity" are an integral part 


of the formulation. 
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The progenitor of this method, as developed by the pres- 
ent investigator, will reproduce all the standard results of 
structural mechanics (with some interesting addtional features) 
when the simplest possible displacement assumptions of the 
hierachy are introduced into a non-classical virtual work 
theorem. The results of this research should be considered 
as a serious candidate to provide the basic equations to be 
used in the next generation of finite element programs for 
routine structural analysis. It may be desireable to give 
these finite element method (FEM) programs self-diagnostic 
capability such that the programs evaluate the input and 
automatically select the proper hierachy of equations that 
will most properly and most economically solve the given 
structural problem. In other words, the self -diagnostic 
capability should be designed to prevent "overkill'' or 
"underkill" in the solution of any given problem. 
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ULTRASONIC NON-DESTRUCTIVE TESTING DEVELOPMENTS 

Senior Investigators: H. F. Tiersten 

P. K. Das 

Since the last report was submitted, trapped energy mode 
focusing transducers have been fabricated by plating a number 
of concentric rings on a PZT-7A plate. The concentric ring 
structure acts as a Fresnel lens and focuses the ultrasonic 
waves. Two transducers, operating in the frequency range 
centered at 2.5MHz, have been fabricated. The first order 
diffraction pattern from a laser probe is being used to study 
the radiation pattern in water at various distances from the 
transducer. From the near field radiation pattern it has 
been confirmed that the acoustic isolation between the rings 
is excellent. 

A doubly focused acoustic imaging system, which incor- 
porates focusing transducers at both the transmitting and 
receiving ends, is being used for imaging of composite mate- 
rials. Scanning is performed using a Z-80 microprocessing 
system which also stores and displays the image. Images of 
the ultrasonic transmission through composite materials, 
which were prepared for tensile testing, are being obtained 
(See Part I, CAPCcMP, section 4). Preliminary results com- 
pare well with those obtained using other testing procedures 
and reveal certain features which have never been observed by 
means of non-destructive testing procedures. 
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The analysis of mode coupling in thickness-extensional 
trapped energy mode transducers with strip electrodes has 
been completed. The frequency spectra occur in closed loops 
due to the large coupling in PZT-7A. The analysis uncovered 
modes with numerous peaks and valleys across the electrode 
in the strongly coupled regions. The lowest complex disper- 
sion curves for PZT-7A plates have been determined and will 
be* employed in th6 analysis of energy trapping in the lower 
frequency range, which has been observed experimentally. An 
analysis of PZT-7A plates driven into thickness-extensional 
trapped energy vibrations by the application of a voltage 
across strip electrodes is being performed. The plate waves 
in the unelectroded region that carry energy away from the 
trapped mode are being included in the treatment. All pre- 
vious treatments were limited to the consideration of an 
eigenvalue problem; the current effort will enable consider- 
ation of the inhomogeneous forced vibration problem. 
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PHYSICAL PROPERTIES OF EPOXY RESINS AND COMPOSITES 

Senior Investigator: S. S. Sternstein 

This project is concerned with the mechanical charac- 
terization of epoxy neat resins, including particularly the 
effects of moisture, temperature and stress histories and 
ultimately the dependence of the properties of high perform- 
ance composites on neat res-in behavior. 

The first phase of this work has involved the applica- 
tion of inhomogeneous swelling theory developed previously 
for filled elastomer systems to the water swelling of a neat 
resin containing isolated carbon fibers. A principal problem 
area has been the development of a suitable consitutive equa- 
tion for describing the thermodynamic interactions of water 
with neat resin and the subsequent mechanical properties of 
the swollen resin. Extensive data from the literature have 
been used to develop a model constitutive equation. Such an 
equation has now been formulated, and it appears to fit the 
available data. Briefly, the model employs a standard lat- 
tice model calculation for the combinatorial entropy of mix- 
ing and a concentration-dependent Flory-Huggins parameter 
for the excess free energy of mixing. 

The mechanical work function required for the inhomoge- 
neous swelling field problem has been developed from the 
Gaussian theory of rubber elasticity and utilizes a concen- 
tration-dependent crosslink density. In effect, this allows 
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the modulus to depend explicitly on water concentration. 

This is to be contrasted with the usual treatment of rubber 
swelling in which modulus and crosslink density do not de- 
pend explicitly on concentration but rather implicitly 
through the volumetric swelling ratio (a simple dilution 
calculation) . It is not suggested that the revised Gaussian 
model has any particular molecular significance; however, it 
fits the data well. One can conjecture that the model is 
that of a two-crosslink density network, one being due to 
the true chemical crosslinks in the epoxy resin and the other 
a physical crosslink network, e.g. , due to hydrogen bonds. 

It is the latter network which is destroyed by water absorp- 
tion and results in reduced modulus. 

The inhomogeneous swelling theory has been revised for 
application to cylindrical inclusions, e.g., carbon fibers 
and sample calculations have been performed using the new 
constitutive equations. It is found that the overall swell- 
ing ratio of the epoxy is little affected by the presence of 
a carbon fiber. However, the local stresses and strains in 
the vicinity of the interface (between the fiber and epoxy 
matrix) are strongly influenced. Thus, a global view of 
swelling in composites, e.g., by weight gain studies, would 
appear as a simple reduction of water absorption in propor- 
tion to the fiber content. However, a local view leads to a 
much deeper insight into the role which water absorption 
plays. For example, it is found that the radial stress at 
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th3 fiber interface generated by inhomogeneous swelling can 
be as high as 50% of the epoxy tensile strength. 

The calculations have also shown that under appropriate 
conditions determined by the constitutive parameters the 
water concentration at the interface is higher than at further 
distance into the matrix. This is in sharp contrast to swell- 
ing in a filled rubber where the swelling at the interface is 
markedly reduced. 

Additional calculations on swelling instability (e.g., 
void formation) are in progress. A publication on the swell- 
ing field calculations is being prepared. 

The second phase of this work involves the dynamic me- 
chanical characterisation of swollen neat resins. Briefly, 
the water -equilibrated samples are being s ijected to rapid 
temperature scans and the. dynamic mechanical loss factor is 
being obtained at frequencies from 1 Hz to 100 Hz. These 
data will be required for further analysis of the water swell- 
ing of epoxy resins. 
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FATIGUE IN COMPOSITE MATERIALS 

Senior Investigators: E. Krerapl 


During this reporting period activities were mostly de- 
voted to the design of a biaxial fatigue test specimen and 

associated fixtures. The Cyclic Strain Laboratory is equipped 
with an MTS tension-torsion hydraulic test system, and we 

plan to utilize this testing machine in our experiments. 

Much thought has been given to means by which we can insure 
that specimens will not fail at the mounting grips when test- 
ed with this equipment. 

The proposed specimen is shown in Figure 31. It is made 
by wrapping graphite-epoxy tape at ±45® around an aluminum 
mandrel. At present four layers of tape are used. At each 
end, 90® layers are built up and subsequently machined to 
the taper shown in Figure 31. 

The specimen is made to fit into the fixture shown in 
Figure 32. The fixture includes a circular main part (1) 
which oolts to the testing machine crosshead. A circular 
cover plate (2) wedges the specimen into the main part through 
the circular plug (3) with square ends which fits into a cor- 
responding opening in the main part (1) and by the split 
wedge (5) and the ring wedge (4) . 

The specimen is mounted by tightening the bolts on the 
circular cover plate (2) . Then hydraulic pressure is ad- 
mitted through a separate part vrhich pressurizes the plug 
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and keeps the wedge assembly under constant pressure. It 
prevents the loosening of the fit due to relaxation in the 
specimen. 

The axial force is transmitted through the wedges. 
Torque transmission relies entirely on friction between the 
plug (13) and the specimen and the split wedge (5) and the 
specimen. 

The plug (3) has a central bore which can be drilled 
through the square ends. This bore will permit the addition 
of internal pressure as an independent load system in addi- 
tion to the axial loading and the torque presently planned. 

Prototype test specimens have been made with no diffi- 
culty, and machining of the fixture is nearing completion. 
Preliminary tests are planned with the MTS tension-torsion 
system in the next reporting period. 

In parallel with the design and fabrication work de- 
scribed above, uniaxial fatigue tests on {0/±452/90}2s and 
{90/+452/0}2 s graphite epoxy 12-ply laminates as reported 
during the last period continued, with special consideration 
of frequency effects. 
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TRANSVERSE THERMAL EXPANSION OF CARBON FIBER/EPOXY COMPOSITES 

Senior Investigators: R. J. Diefendorf 

C. LeMaistre 


It is important to be able to predict physical proper- 
ties of laminates from constituent properties by means of 
micromechanical analysis. Several authors have developed 
micromechanical theories and have been able to obtain reason- 
able correlation with experimental data. This is particular- 
ly true of elastic and thermal properties with isotropic con- 
(1-5)* 

stituents' . Moreover ^ the longitudinal TCE, de- 

rived by Greszczuk^^^ and Schapery^^^ is even effective for 
composite systems which have anisotropic fibers .graphite/ 
epoxy) . Their expression for is: 


. _ “mVm * 

- Vj) + Vf 


( 1 ) 


where a, V, E stand for TCE, volume ratio and Young's modulus; 
and m and f subscripts stand for matrix and fiber, respective- 
ly. This equation assumes that Poisson's ratio of the fiber 
and matrix are equal. The good agreement is obtained because 
the fiber properties dominate in the 0* orientation, and 
fiber properties are measured parallel to the fiber axis. 

Prediction formulas for transverse TCE do not fare as 
well for graphite reinforced composites as they do for other 
composites due to the anisotropy of the fibers. The theories 


* 

Nos. in raised parentheses refer to the References liisted on 
pages 98 and 99. 
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to be presented in the following paragraphs make the follow^ 
ing basic assumptions: 

1) The ply is macroscopically homogeneous, linearly 
elastic and generally orthotropic. 

2) The fibers and matrix are linearly elastic, homo- 
geneous and isotropic. 

3) There is complete bonding at the interface of the 
constituents, and there is no transition region 
between them. 

4) The ply is initially in a stress-free state. 

f 8) 

5) The fibers are regularly spaced and aligned. 

A mechanics of materials approach was used by Greszczuk^^^ 
to determine and ^ 22 * quantity was presented in 

Equation 1 and 6^22 expressed as follows: 

E e + a E (1 - 3)] (2) 

zz E 2 00 mm 

a = a_(l - 23 ) + 2a-3 - v (a- - a ) (1 - 23) (3) 

o m r m r ra 

E E^ 

p = m_f 

^o tE.d - 23) + 2E^3l 
r m 

E» = E 3 + E (1-3) 

2 o m 

../3 

(7) 

This formula assumes circular filaments. Schapery devel- 
oped bounds on the TCE of isotropic constituent composites by 
employing extremum principles of thermoelasticity. As shown 
before. Equation 1, Schapery *s formula for is the same as 
that derived by Greszczuk. Schapery also presents a simple 


(4) 

(5) 

( 6 ) 
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expression as an upper bound for ^22 derived by Wang Fo- 

( 5 ). 


Fe 




The results from formula (7) compare well with experimental 

results from fiber glass composites. 

Pahmy and Ragia-Ellozy and Kahn^^^^ applied finite 

element analysis to determine and for a lamina. This 

11 22 

approach provides a rigorous solution for ^ 22 ^ but requires a 
sophisticated microstructural analysis of the composite sys- 
tem and does not provide a readily accessible equation for 
determining ^22 from constituent properties. Kahn also de- 
velops an equation for ^ 22 ' treating the fiber as a solid rod 
of radius 'a' and the matrix as a tube of inner radius 'a' 
and outer radius 'b*. The matrix tube is subjected to an in- 
ternal pressure, P, whose sign depends on the thermal mis- 
match between fiber and matrix: 


®22 ' “m 


AT 


- V a,, 1 

0m D E 1 - X zm 

m 


( 8 ) 

(9) 


E,(l - V ) 
r m 


V. 


'^z P[1 + E (1 - V ) (1 + V.) + E V (1 + V )^1 - V. 

m mf fffm f 

These equations result in a curve which falls between the rule 

of mixture lower bound solr.tion and Schapery's upper bound 

equation. Finally, Levin provides an anisotropic analysis 

for the thermal expansion coefficients. No comparison is made 

to experimental results for almost all analyses. 
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In an attempt to evaluate the applicability of these 
theories, thermal expansion specimens were fabricated and 
tested on a push-rod dilatometer. The material selected 
for this investigation was the Piberite 1048 AE graphite/ 
epoxy system. Components consist of Union Carbide T-300 
graphite fiber and Piberite 's 121*C (250®P) 48A resin in 
prepreg form. Test specimens were fabricated on a Dake 75- 
ton press with 19'! x 19" programmable heated platens. The 
cure cycles were programmed with a Data Trac and fully moni- 
tored. Nine specimens, 6" x 6", were fabricated with vary- 
ing fiber volume, from which five expansion samples, 2" long, 
approximately 1/2" in highth and thickness varying depending 
on fiber volume, were cut. All nine specimens contained 32 
plies unidirectional prepreg cured at various pressures, 
quantity of bleeders and excess resin to achieve the required 
resin content. Neat resin was cured in the same manner as 
the laminates under vacuum pressure to remove porosity. Sim- 
ilar dimensions were used for the resin TCE specimens. Dog- 
bone tensile specimens were fabricated for determination of 
tensile modulus and Poisson's ratio of the resin. 

TCE specimens were conditioned in an air circulating 
oven at 110®C for 24 hours and stored in a desiccator until 
tested in an attempt to eliminate moisture induced variations 
in the thermal expansion measurement Specimens were 
v;eighed before and after testing to determine moisture loss 
during testing. 
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Thermal expansion measurements were performed on a modi- 
fied Thermo-Physics Corporation push rod dilatometer. Sam- 
ples were heated from room temperature to 100*C at a low rate 
of 0.75*C/min. to minimize discrepancies which may arise from 
the variation in thickness and thermal conductivity of the 
different fiber volume samples. The linear TCE was derived 
from a continuous recording of length, utilizing a precision 
LVDT, with resulting accuracy of 0.1 x 10”V®C. 

Two methods were used to determine fiber volume; 1) vis- 
ual count method using exploded photo-micrographs of the com- 
posite microstructure and 2) Boeing's MS8-212A nitric acid 
digestion method. Densities were determined utilizing ASTM 
D792 "Tests for specific gravity and density of plastics by 
displacement. Method A". 

Mechanical testing was performed on an Instron screw- 

type universal testing machine at a crosshead rate of 0.005"/ 

min., and strain data was recorded on a Data Logger. The 

1048 AE composite's mechanical property data used were re- 

(14) 

suits of early testing by Kowalski 

The results of individual testing indicate that the 24 
hour thermal conditioning of the TCE specimens was effective 
in minimizing the moisture response on expansion. The typi- 
cal peaking, probably caused by the T depression, evident 
(X2 13 15) 

in past work ' ' ’ due to absorbed moisture, was not 

noted. Weight loss related to moisture desorption during 
testing was less than .03% in all samples tested. 
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Table XII shows average transverse TCE from 30*C to lOO^C 
versus fiber volume. These results are plotted on Figure 33. 
From these data points, an estimate of ^22 individual 

fiber was calculated by utilizing a least-mean-squares extrap- 
olation to Vj “ 100%. It is noted that a modified rule of 
mixtures behavior is assumed and that a straight line is valid 
to approximate ^ 22 * From this analysis it is estimated that 
^22 T-300 graphite fiber has a value of 8.8 x 10”V°C in a 

temperature range of 30®C - 100®C. This result is similar to 
. (15) 

a result obtained by Ishikawa who estimates a value of 8.6 
X 10 ^/®C. However, it is important to note that Ishikawa 's 
method appears to be self-fulfilling in that 6(22 is determined 
by extrapolating values from a theoretical curve and plugging 
the results back into the theory to prove its fit. 

In order to correlate this estimated value of 6^22 of T- 
300 to Schapery s upper bound solution, elastic constants of 
the matrix are required. Test results indicate that = 

g 

.55 X 10 psi and = .35. These values agree with typical 

( 16 ) 

values of and for other epoxy resins used in composites , 
The upper bound equation presented by Schapery was modi- 
fied 


a 


22 


+ (1 + V 

m mm 


fl 2 )“f 22 ^f ‘ ^m'^m^ (7a) 


to include fiber anisotropy, a. is defined by Equation (1) 

Oi 

and would appear to be close to the value a^ for high modulus 

^11 

fibers at moderate or high volume fractions. Direct applica- 
tion of this equation gives a theoretical curve that is higher 
than all experimental values, ; because the boundary condition 
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TABLE XII 


Specimen 

Number 

TEST RESULTS : TRANSVERSE 

TCE VERSUS FIBER VOLUME 

Density 

(g/cc) 

(nitric digest) 

(micro-photo) 

TCE (xlc'^^C) 
(ave. : 30-100*0 

11 

1.542 

67.2 


23.5 

12 

1.595 

70.1 


27.2 

13 

1.454 

— 


25.1 

14 

1.586 

64.2 


26.6 

15 

1.585 

69.2 


25.2 

ave. 


67.7 

68.7 

25.5 


21 

1.553 

— 


30.4 

22 

1.563 

65.1 


29.7 

23 

1.524 

64.7 


29.7 

24 

1.533 

64.6 


28.2 

25 

1.580 

69.3 


26.8 

ave. 


65.9 

63.0 

28.9 


31 

1.628 

67.3 


31.9 

32 

1.550 

63.9 


28.7 

33 

1.543 

66.1 


30.1 

34 

1.514 

62.4 


33.7 

35 

1.542 

63.8 


31.0 

ave. 


64.7 

64.2 

31.1 


41 

1.530 

61.9 


32.3 

42 

1.557 

62.1 


32.9 

43 

1.579 

62.9 


29.7 

44 

1.557 

62.4 


33.9 

45 

1 . 53C 

63.8 


31.3 

ave. 


62.3 

60.9 

32.0 


51 

1.532 

59.2 


36.2 

52 

1.537 

59.4 


31.7 

53 

1.525 

58.8 


35.4 

54 

1.500 

57.4 


34.7 

55 

1.513 

58.1 


33.4 

ave. 


58.6 

56.3 

34.3 
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TABLE XII continued 


Specimen 

Number 

Density 

(g/cc) 

(nitric digest) 

(micro-photo) 

TCE (xl0”V®C) 
(ave.:30-100‘»C) 

G1 

1.470 

43.0 


34.7 

62 

1.451 

50.9 


37.1 

63 

1.443 

50.4 


39.9 

64 

1.454 

51.2 


40.5 

65 

1.449 

51.3 


39.1 

ave. 


49.4 

51.5 

38.2 


71 

1.501 

52.1 


38.3 

72 

1.492 

54.4 


36.6 

73 

1.489 

54.4 


38.1 

74 

1.489 

54.1 


39.1 

75 

1.478 

54.9 


35.0 

ave. 


54.9 

53.6 

37.4 


81 

1.354 

46.9 


40.0 

82 

1.401 

46.8 


42.4 

83 

1.407 

— 


38.6 

84 

1.417 

48.1 


41.2 

85 

1.431 

52.0 


38.8 

ave. 


48.4 

49.6 

40.2 


91 

1.393 

40.4 


46.1 

92 

1.343 

38.4 


44.1 

93 

1.328 

37.3 


41.0 

94 

1.342 

38.7 


44.4 

95 

1.389 

41.3 


43.3 

ave. 


39.2 

36.6 

43.3 


101 

1.185 


51.9 

102 

1.201 


53.9 

103 

1.200 

00.0 

52.5 

104 

1, 2'! 3 


52.8 

105 

1 ■>0L 


54.8 


ave 


53.2 


TCE versus Fiber Volume 
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at » 1 is not satisfied. In order to satisfy this condi> 
tion, the following must be true: 


lim d,, « 

Vj-^1 

(1 -r V- )a- - d 

ti2 122 ® ^12 

(11) 

-I, d22 

a- - V- - d V- 

^22 ^22 ^37 ^22 ’ ^12 

(12) 


*> d. = a,' 

U ^ 0 0 

(13) 


The resulting equation is: 

a 


‘22 


+ (V- - V V )a. 
m m m r m m r 


(7b) 


22 


Input values used in the theoretical prediction plotted on 
Figure 33 ares: 


v„ = .35 
in 


a„ = 53.2 X 10“V*C 
m 


E = .55 X 10” psi 
m 


^12 


22 


= 8.8 X 10 " 5 /»C 


= 32.5 X 10 psi 


It is seen from Figure 33 that excellent correlation exists 
between experimental data and the results calculated from 
Equation (7b) . The reason for the need to modify the form 
of the equation falls back on the hidden assumption in the 
derivation of Equation (7a) in the third term that the con- 
stituents are isotropic in nature. Equation (7a) assumes 

that 61 = 6 i, , from isotropy and d " where K is the 

a 11 a K 

volume ratio of the bulk modulus (i.e., EK V^) and K is de- 

n n 


finod to be: 


^11 ^22 ^33 " 


(14) 


Ka is defined in the same manner. In Equation (14) , is 
equal to + 6122 ^ 33 ' in the case of isotropic 
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constituents can be reduced to Substituting this 

isotropic constituent condition in Equation (14) and assuming 

that — ~ Equation i7) is derived by energy principles. 

E K 
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